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Merging the benefits of an NLO calculation with event generators is a topic of current interest. 
Among other things, it is advantageous to be able to generate unweighted events according 
to NLO matrix elements. I report on an algorithm which generates a sequence of unweighted 
momentum configurations, each configuration containing either n or n+1 four-vectors, such 
that for any infrared safe observable the average over these configurations coincides with the 
NLO calculation up to errors of a resolution variable. 

1 Introduction 

As high energy experiments move towards a new generation of colliders, accurate simulation of 
collision events becomes increasingly important. In particular this is the case for QCD, which 
often constitutes the bulk of the events. Theoretical predictions are in most cases either based 
on exact perturbative calculations, resummations or event generators. These methods have 
complementary regions of applicability. Exact perturbative calculations give a prediction for an 
observable at the parton level as a power series in the coupling constant. The state-of-the-art 
are numerical next-to- leading order programs, which allow the calculation of any infra-red safe 
observable for a given process. They give a reliable prediction, if all scales entering an observable 
are of the same order of magnitude, but are limited to a few partons only. If an observable 
depends on a large ratio of two scales, large logarithms appear in the perturbative expansion. In 
this situation one often reorganises the perturbative expansion and resums the large logarithms 
to all orders. Within resummation, the state of the art is next-to-leading logarithmic accuracy. 
Up to now, these resummations have to be carried out for each observable individually. Event 
generators follow a different philosophy and rely on a three-stage process: Starting from a hard- 
scattering process, parton showers are generated in the second stage. Finally, a hadronization 
model converts partons into hadrons. The first two stages, hard scattering and parton showering, 



can be described in perturbation theory, while hadronization involves long-distance physics 
and is non-perturbative. Event generators usually only employ leading-order matrix elements 
and partons showers accurate to leading logarithmic accuracy. They describe reasonably well 
multiple collinear emission and they generate unweighted events, since they only use leading- 
order matrix elements. However, they do a rather poor job for multi-jet events. This is due 
to the fact that they only use the simplest 2—5-2 and 2^3 LO matrix elements, that the 
showering algorithm does not populate the complete phase space and that color interference is 
only approximated by angular ordering. Since multi-jet events become increasingly important for 
new physics searches, it is a topic of current interest to merge the benefits of exact perturbative 
calculations and event generators. Several issues have to be addressed for an event generator 
according to NLO matrix elements: 

• Double counting: Gluon emission can originate from the hard scattering as well as from 
parton showering. 

• Color flow: Matrix elements are usually color summed, where as the algorithms for show- 
ering and hadronization are based on open color lines. 

• Negative weights: A standard NLO program generates n- and (n + l)-parton phase space 
contributions, which are finite, but not necessarily positive definite. 

Recent developments in this field considered combining parton showers with multi-parton tree- 
level matrix elements 1-3 , with resummations of soft gluonsu, and with NLO matrix elements 5-9 . 
In this talk I focus on the issue of negative weights. Negative weights are not a problem per se, 
but can seriously waste computer resources, if events are passed through a detector simulation. 
It is generally the case that the CPU time needed for detector simulation exceeds the one needed 
for the generation of an event. It is therefore desirable to avoid them from the very beginning. In 
standard NLO calculations negative weights occur since individual contributions to an observable 
are not necessarily positive definite. It is therefore natural, to consider first a "partonic event 
generator", leaving parton showers and hadronization for the moment aside. Here I report 
on an algorithm which generates a sequence of unweighted momentum configurations, each 
configuration containing either n or n+1 four-vectors, such that for any infrared safe observable 
the average over these configurations coincides with the NLO calculation up to errors of a 
resolution variable £9. I focus on electron-positron annihilation with massless partons. The 
algorithm is independent of the specific hard process and can be implemented on top of existing 
NLO programs. After an outline of the basic idea, I focus on a few technical challenges which 
have to be addressed in order for the algorithm to run efficiently. 

2 The algorithm 

Before stating the algorithm I briefly recall the main features of standard NLO calculations. 
The calculation of an observable at NLO receives contributions from the Born matrix element, 
the virtual corrections and the real emissions. Taken separately, the last two contributions are 
divergent, only their sum is finite. The subtraction- or phase space slicing method can be used 
to render them finite individually. Within the subtraction method one adds and subtracts a 
suitable piece du^llil: 




The modified virtual corrections and the modified real emissions are now separately finite, but 
not necessarily positive definite. The basic idea to generate unweighted events according to NLO 



matrix elements is relatively simple and inspired by a slicing approach (although it will entirely 
based on the subtraction method): Consider a (n + l)-parton configuration. If all partons are 
separated by at least y res according to a resolution criteria, one returns this (n + l)-parton 
configuration. If on the other hand a parton is not resolved at a scale y re s, one averages over the 
real emissions up to the scale y re s, combines the result with the virtual corrections and returns 
a n-parton configuration. By choosing y res large enough the sum will have a positive weight. 
For an efficient implementation I would like to discuss three technical points: 

• A modified version of the subtraction method, such that the subtraction terms are non- 
vanishing only for y < y res - 

• A parameterization of the phase space, which starts from a n parton configuration and 
inserts an additional soft or collinear particle, such that the integration over the unresolved 
part involves only a three-dimensional integration and, in addition, such that the dominant 
1/sjj-behaviour of the integrand is absorbed into the measure. This ensures that the 
unresolved region can be integrated efficiently. 

• The Metropolis algorithm to generate events with constant weight according to a (com- 
plicated) probability distribution. 

I start with a modified version of the subtraction method: The subtraction method has in general 
the advantage that it does not introduce additional approximations to the matrix elements or 
to the phase space. The only approximation made for the algorithm presented here, consists in 
replacing a (n + l)-parton configuration by a n-parton configuration, whenever two partons are 
not resolved within y res - In order for the basic idea outlined above to work, nothing should be 
subtracted if all partons are separated by y res - In the original version of the dipole formalism 
1111 the approximation term is given by a sum over dipoles J2^ij,k and subtracted over the 
complete phase space. However, it turns out that the approximation term can still be integrated 
analytically if the region where it is subtracted, is restricted 111 : 

d<J A = ^2 T> ij:k Q (y ijik < y res ) (2) 

Here, = 2piPj/ (2piPj + 2pip k + 2pjPk) is the appropriate resolution variable if partons i and 
j become close to each other and k is the spectator. 

The second point is a suitable parametrization of the phase space. For any given "hard" n-parton 
configuration one needs a parametrization for the insertion of an additional parton, which can 
become soft or collinear to one of the n hard partons. This is done as follows £3: The phase 
space measure of the (n+ l)-parton configuration can be related to a n-parton configuration by 

ji , , ds as ds sb d(j) s 

= dK 4(2vr)3, afe • (3) 

The invariants s as , s s b and the angle 4> s are parametrized in terms of three random variables v\ , 
V2 and t>3, uniformly distributed in [0, 1]: 

( s min\ 1 ( s min\ 2 , n /,\ 
, S s b = Sab ,9s = 2ttv 3 . (4) 
Sab ) V S a b J 

Smin is an internal technical parameter, which is chosen to be small. It is possible to recon- 
struct the four- vectors of the n + 1 parton configuration from these invariants. The phase space 
measure becomes 

dcf> n+ i = d^n -j- „ SasSsb l n 2 (fUEL) dvidv 2 dv 3 (5) 
4(2tt) 2 s ab V s ab J 



and the 1/sy-behaviour of the integrand has been absorbed into the measure. 
With an appropriate choice of y res one obtains now events with positive weights. In general, 
the probability distribution will be a rather complicated function and unweighting according 
to the standard acceptance-rejection method might not be feasible. The method of choice is 
therefore the Metropolis algorithm, which allows to generate unweighted events according to a 
multi-dimensional probability density P{x\, ...,Xd), which not necesarrily factorizes. One starts 
from a state x = (x\, ...,Xd) and randomly generates a new candidate x! . One then calculates the 
quantityAS = — \n(P{x')/P{x)) and accepts the new candidate with probability min(l,e~ AS ). 
The way new candidates are suggested is arbitrary, restricted only by ergodicity and detailed 
balance: Each state must be reachable within a finite number of steps and the probability of 
suggesting the state x^ given we are in state x\ has to be equal to the probability of suggesting 
x\ given we are in x^. 

To summarize, the complete algorithm to generate unweighted events according to NLO matrix 
elements is therefore as follows: One first picks a set of random numbers U\, ...,u m ,x, 1)1,1)2,1)3, 
where u\,...,u m define the "hard" n-parton configuration, x decides where the "soft" particle 
gets inserted and v\ , x>2 , V3 (together with the "hard" configuration and the information where 
the "soft" particle gets inserted) yield the "soft" {n + 1) event. One then checks for the (n+ 1)- 
parton event if all are greater than y res , e.g. if all partons are well separated. If this is 
the case the probability P' is proportional to the real emission matrix element and the (n + 1) 
parton configuration will be returned, if the event is accepted by the Metropolis algorithm. If 
not all partons are well separated, the "hard" n-parton configuration has the probability 

P' ~ \M^J 2 + 2ReM^ n *MUoo P + (M n Born \I\M n Born ) 

~ V ij,k® (Vij,k < yres)) (6) 

The last terms in the first and second line are the added and subtracted dipole terms, which 
ensure that the expressions over the n- and (n + l)-parton phase space are finite separately. 
The second line involves a three-dimensional integration over the unresolved phase space. In 
both cases, e.g. if the event is classified as a n- or (n + l)-parton event, the new candidate 
gets accepted with Metropolis probability min(l, P' /P), where P' is the probability of the new 
candidate and P is the probability of the previous event. 

3 Outlook 

I presented an algorithm to generate unweighted partonic events according to NLO matrix 
elements in electron-positron annihilation. An important ingredient is an efficient integration 
over the unresolved region on a point-by-point basis in hard phase space. The extension to 
initial-state partons or massive partons should be feasible. 
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